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ABSTRACT. In a previous paper a technique was developed for finding reconstruction 
algorithms for arbitrary ray-sampling schemes. The resulting algorithms use a 
general linear operator, the kernel of which depends on the details of the scanning 
geometry. Here this method is applied to the problem of reconstructing density 
distributions from arbitrary fan-beam data. The general fan-beam method is then 
specialized to a number of scanning geometries of practical importance. Included 
are two cases where the kernel of the general linear operator can be factored and 
rewritten as a function of the difference of coordinates only and the superposition 
integral consequently simplifies into a convolution integral. Algorithms for these 
special cases of the fan-beam problem have been developed previously by others. 

In the general case, however, Fourier transforms and convolutions do not apply, and 
linear space-variant operators must be used. As a demonstration, details of a fan- 
beam method for data obtained with uniform ray-sampling density are developed. 








REVIEW. 


In a previous paper [1], I developed a technique for finding reconstruc¬ 
tion algorithms applicable to arbitrary ray-sampling schemes. This general 
method was applied to the problem of reconstruction from parallel-beam data 
with uneven spacing between rays and uneven spacing between projections. The 
results were based on Radon's famous integral [2], 

f (r,<{>) = — j J (- £) ~ p(*,e) d£ de (1) 

T o 

where pU,e) is the density integral or ray-sum measured along the ray inclined 
e with respect to a vertical axis and passing within a distance i from the cen¬ 
ter of the region being scanned (see figure 1). Further, f(r,4>) is the density 
at the point with polar coordinates (r,<j>) in this region, while t = s. - r cos(e 
is the perpendicular distance between the ray and this point. 

When ray sums in a given projection are spaced evenly in a, and projections 
are spaced evenly in e, a simple reconstruction method can be found directly 
from equation 1 by approximating both integrals by sums and approximating the 
partial derivative by an appropriate first difference. 



t 
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UNIFORM SCANNING COORDINATES. 


When spacing is uneven, it is helpful to introduce first new ray-sampling 
coordinates 5 and n> chosen so that successive rays in a generalized projection 
correspond to evenly spaced values of 5 , while successive projections corres¬ 
pond to evenly spaced values of n- Radon's integral can then be transformed 
to this new coordinate system using the Jacobian, 


j _ 3 £ _ 39 


32. . 30 


35 3n 3n 3C 


and equation 1 becomes. 


f(r,<{>) = f J (" J(C»n) d5 dn 

4ir 2 




It is possible to show that this can be rewritten as 


f(r,<f>) 



3D 30 3D 30 -i ,, r . 

. - - —C. . - f|£ df, 

35 3n 3n 35 



It is not clear whether this forms a good basis for a reconstruction algorithm 
in the general case, since it seems to imply that computations must be carried 
out across projections as well as within projections. 
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GENERAL PARALLEL BEAM METHOD. 

In the previous paper [1], the emphasis was on parallel-ray scanning, and, 
in this case, £ is a function of z only, while e is a function of n only. The 
Jacobian then reduces to. 




and equation 4 simplifies as follows, 

^ / c / ( -H? (6 > 

Here t = z - £’, where z - n(s), while z' = z(z'), and z' is the value of z 
associated with the ray that passes through the point (r,<f>). The above can 
be conveniently split into an outer and an inner integral: 


f(r,<*>) = —— ( gU‘>n) dn 

4tt2 J 9n 

g(5 >n) J £( 5 ) - z(z’) I? ^ 

If these integrals are approximated by sums, one obtains: 

f(r,<{>) : 4 ^- I (C) 50^ 


(7) 

( 8 ) 
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q., . 

J 


= - I 


(p ij " P (i - 1)j ) 


(IT 


^TTT) 


(io) 


This straightforward set of equations is one result of the analysis in the 
previous paper (equations 29 and 40 in [1]). Here fie. = (e^ + -j “ 0j _ -j 


)/2 


th 


is the angular interval associated with the j ul projection, while fci - 


„.th 


(s,. + i. , , )/2 is the value of a corresponding to the center of the i 1 "" beam 
i i + l 

XL 

The left edge of the beam striking the i detector corresponds to and the 
right edge to £... + -j (see figure 2). The density integral obtained from the 
i** 1 detector in the j^* 1 projection is p.... 

Finally, note that g.(z') has to be found by interpolation from the dis- 
crete set of values, {g., -}. If linear interpolation is to be used, one can 

J 

work with the values g^ and g^., + -|)j» where 


1 lU') <^1 + 1 


( 11 ) 





RELATION TO CONVOLUTION-BACKPROJECTION METHOD. 

By splitting the second sum and rearranging its terms, one arrives at an 
alternate form (equation 38 in [1]), 




_4. 

^i’ + 1 




That is, the sequence {g..} is obtained from the sequence fp^^> by a general 
linear operator. This is similar to a convolution except that the weights or 

filter coefficients are spatially variant. 

One has only to fix the width of the detectors, at SJt say, to be able to 
relate this result to the well-known convolutional-backprojection method. In 

this case, 




This amounts to convolution of {p- -> with a filter function F^, where 

»J 


p =-^- for k f 0 and F = 4. 

k 4k 2 - 1 0 



This happens to be the particular set of filter coefficients popularized by 
Shepp and Logan [3]. 

My previous paper contains other formulations for this problem as well as 
a simulation of the method for reconstruction from ray sums collected with 







uneven spacing [1]. The main point is that convolutions are inappropriate in 
the general case and must be replaced by spatially varying operators or super 
position integrals. 
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RECONSTRUCTION FOR ARBITRARY RAY-SAMPLING SCHEMES 


In the general case, equation (4) does not seem to provide a good starting 
point for analysis. Instead it is helpful first to remove the partial deriva¬ 
tive from Radon's integral by partial integration. This has to be done care¬ 
fully since the inner integral is singular. It can be written as 


lim 

e-*0 


l -£ 


- I) fx pU,e)d« + / 


(- PU,9) d£ (15) 


& * +£ 


Integrating by parts one obtains (equation 8 in [1]), 


C JT 00 

= 4 ^ / ^0 / F e (t)p(fc,e)d£ de 


where 


F (t) = 1/e 2 for |t| < e 
0 


(17a) 


= - 1/t 2 for |t| * e 


(17b) 


Introduction of the transformation to uniform scanning cbohdinates (c,n) now 


leads to 
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f(r,$) = — 

4*TT 2 

g(r,<f>,n) = f F (t)J(5.n)pU,n)de O 9 ) 

This forms the basis for reconstruction methods for arbitrary ray¬ 
sampling schemes. In this paper the concern will be with techniques for re¬ 
constructing density distributions from ray sums collected using fan beams 
(see figure 3). Modern apparatus for computerized tomographic analysis 
typically produces projection data in this form and there is a practical need 
for accurate and rapid reconstruction methods for a variety of different 
schemes for sampling the fan beam. Such methods had been found previously for 
two special geometries [4, 5, 6]. Here techniques will be developed that cSh 
be used for arbitrary ray-sum collection schemes. 


I 


g(r,$»n) dn 


(18) 
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SOURCE POSITIONS DISTRIBUTED AROUND THE CIRCUMFERENCE OF A CIRCLE. 

Let the source be located at (D,tt/ 2 + g), where D is the radius of the 
circle (see figure 1). Let a ray be emitted in a direction that makes an angle 
a with the source-to-origin line. Clearly a and g are as good for specifying 
a particular ray as 1 and e are. For fan beams, these new parameters will be 
more directly useful, and so the relationships between the two sets of variables 
will be needed. From figure 1, 

a = D sin a and 9 = a + g (20) 

a = sin" 1 (£/D) and g = e - sin _1 (£/D) (21) 

If £ and n are uniform scanning coordinates, then it is natural to let 

a = a(c) and g = p(n) (22) 


where a and g are continuous, differentiable monotonic functions of £ and n. 

^ . * ■ * • 

Then (see equation 2), 


i - 1L 1 ® 
d " 3£ 3n 


since 



(23) 


Further, since 9 = a + g. 


32, 3a 3g 
3a 3£ 3n 


3g 

9n 


D COS a 


3a 

9C 



Now J is a factor in the inner integral (equation 19), but the first term of the 




above product can clearly be brought out of the inner integral and incorporated 
into the outer integral (equation 18). The last term of the product will depend 
on the way in -which rays in a fan are sampled. This corresponds to the place¬ 
ment of detectors in the fan and depends on the scanning scheme used. We will 
study several cases after developing a few more tools that will be needed. 



THE PERPENDICULAR DISTANCE FROM A POINT TO A RAY FOR FAN BEAMS. 

From the diagram (figure 4) we can develop a useful new way of writing t, 
the distance between a given ray (a, 3 ) and a given point (r,<j>): 

t = K sin (a - a') (25) 

where, by the cosine rule for triangles (see figure 5), 

K 2 = r 2 + D 2 + 2r D sin (3 - <J>) (26) 

Here a' is the value of a corresponding to a ray from the source which passes 
directly through the given point (r,<f>). Note that K is simply the distance 
from the source to the point (r,<j>) and thus clearly does not depend on a. 

From the diagram we can calculate a 1 as follows (see figure 5): 

K sin a' = r cos (3 - <J>) and K cos a' = D + r sin (3 - 4>) (27) 


and so. 


r cos (6 - <i>) 

[D + r sin (3 - 4>)3 


tan a' 


(28) 
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SOME PROPERTIES OF THE FILTERING FUNCTION. 

Note that if c f 0, then 

F e (ct) = £7 ^ z / c (t) (29) 

. _ . ' f ( _ . 

* 

This result can be easily checked by separately considering the cases [ct| < e 
and )ct[ * e. From this it follows that 



(30) 


and 


lim f F (ct) a (t) dt = X lim f F (t) a (t) dt (31) 

e ->• 0 J e c e + 0 J e 

— 00 —CO 


Furthermore if |b(t)| k b > 0 is continuous and differentiable with respect to 


t, 

' im f F c [b(t)t] a(t)dt = ,1m f F t (t) f||iy dt (32) 

® ® -oo £ 0 _ oo 

These results are useful in deriving the general reconstruction formula for fan- 
beam scanning schemes. 
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RECONSTRUCTION ALGORITHM FOR ARBITRARY FAN BEAM GEOMETRIES. 

Using the general result for reconstruction (equations 18 and 19) based on 
Radon's inversion formula [2] and using the expressions for t and J just de¬ 
rived, we have 

f ( r »<i>) = 4 ^t j 9'(r,<M) d n 

lim f 

g 1 (r ,(f> ,n) = / F [K sin (a 

e + 0 J c 


(33) 

- a')]D cos a || p(?,n) d£ (34) 


According to a result just developed (equation 31), the factor K can be extracted 
from the argument of the filter function F. Since K does not depend on a it 
can be further removed from the inner integral (equation 34) and placed into 
the outer integral (equation 33). The inner integral then no longer contains 
terms which depend explicitly on r and <j>, only terms which are a function of 
a'. So the above can be rewritten. 


f(r 


,0 = 4^/ 9"(«\ 


n) (1/K 2 ) f^dn 


(35) 


lim f « 

”{a\n) = / F [sin (a - a')] D cos a ff- p (£,n) d£ 

e - 0 J e 35 


(36) 



From the fact that the inner integral is a function of o', and not r and 
4 ) explicitly, we conclude that the reconstruction algorithm can be arranged 
efficiently. That is, for all fan beam schemes with source positions on the 
circumference of a circle and sampling of the fan independent of source position, 
one need not explicitly calculate the contribution of each ray to each point in 






the reconstruction. 

We are now ready to develop specific reconstruction methods for a variety 
of fan-sampling schemes. Some special schemes will lead to simplifications 
which can be interpreted as pre-multiplication, convolution and post-multi plica 
tion. In general, however, the inner integral remains in the form of a general 
linear operator or superposition integral. 
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UNIFORM SAMPLING ALONG A LINE AT RIGHT ANGLES TO THE SOURCE-TO-ORIGIN LINE. 


Rays are sampled evenly in X (see figure 6), so it is natural to let 

C = A = D tan a (37) 

D cos a (||) = cos 3 a (38) 


Also 


sin a = 


/W'+z* 


and cos a = 


D 


SWT1? 


(39) 




D cos a(||) = (40) 

Further, if we let a' be the value of a corresponding to the ray through the 
point (r»<f>) and s' the corresponding value of s> then 


sin a' = 


/ W + V 7 


and cos a 1 = 


D 




(41) 


So, 


sin (a - a') = sin a cos a' - sin a ' cos a (42) 
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sin (a - a‘) = 


D 


/ D 2 + 5 2 / D 2 + l 


rz 


(? - S') 


(43) 


Hence, 


g"U',e) = 


lim 

E 0 



1 U - £') D 


m + v 7 


/D 2 + 5 2 


D 3 


[D 2 + 5233/2 


pU,B)d€ 


(44) 


We can move the multiplier of (5 ~ £') out of the filter function F to get: 


g"(C,&) = (D 2 + V 2 ) 


lim 


e 



D 


F e (« - V) P(«.6) dc 

e /D 2 + 5 2 


(45) 


Note that (from equation 37 or 41), 


5 ' = D tan a' 


(46) 


and so 


(D 2 + c' 2 ) = D 2 sec 2 a' 


(47) 


Also, combining this with the term (l/K 2 ) in the outer integral one gets 


(D 2 + 5 ’ 2 ) D 2 _ D 2 

Kr = IK cos V] 2 = [D + r sin (e - *)] 2 


(48) 


and so, finally. 



"A +<*rm ry*:- j.*.- ^ r 






(49) 


2tt 

f ( r,<f> ) = 4^ f 9 "' ^ X ' * [D + r sin (3 - 4»)]^ dB 

o 

n * 4*oo 

lim r n 

S"'(X , .B)- / F (X - A') — f -- pU.B) dx (50) 

e * o-' 

— 00 

Here finally we have used the coordinate natural to this particular scanning 
scheme, namely x, the distance from the origin at which the ray intersects the 
line drawn perpendicular to the source-to-origin line. 

It is important to note that in this case the argument of the filter function 
F contains only the difference of the two parameters X and x'. The above is 
thus almost like a convolution, except that one has to pre-multiply the projec¬ 
tion data p(x,e) by a factor depending on the position of the ray in the fan. 
Similarly, the convolved data, g"' , is used in the outer integral after post¬ 
multiplication by a factor which depends on the position of the point (r,<j>) i 
the fan currently being processed. 
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DISCRETE APPROXIMATION. 

Finally, we have to approximate these integrals by sums because in practice 
only a finite set of ray sums is available: 


f(r,*) = 4^7 £ 9 j (x,) 7^ 


D* 


[D + r sin (3. - <(>)] 


2 66 i 


(51) 


Q • i • = / i F * * i 

3 1'j . 1-1 

J 1 


D 


P.. 6X 


/ b'2~T~X~2- 


(52) 


Here 


D r cos ( 3 . - <|>) 
A' = D tan a' = pTT r si - nTe 7" - ~iT 

J 


(53) 


_.th 


■ th 


In the above set of reconstruction equations, P.., is the i v " ray sum in the j 
fan, while 6\ is the (fixed) interval between intersection of successive rays 
with the line perpendicular to the source-to-origin line. The angular interval 
associated with a particular fan is 53 ., where 

J 


< 5 f>j = (~ $j_'| )/^ 


(54) 


The filter factors are 


W k 

F k (k 6 x) T 


k j 0 


(55) 


co 


0 - 2 r k 


F_ = - 2 V F 

k = 1 


(56) 







As mentioned in the paper on which this analysis is based, the weights w are 

k 

chosen to provide good numerical approximations for the singular integral. 

Typical choices are: 

1. = 2 for k even, w. = 0 for k odd 

2 . w k = 4k 2 /( 4 k 2 - 1 ) 

3. w = 1 

k 

Note that the above operation can be viewed as a pre-multiplication of the 
ray sums by D/CD 2 + x.j 2 ] 1/2 , followed by a convolution, with a final post-multi¬ 
plication by D 2 /[D + r sin ( 3 . - $)] 2 . While the method is not strictly con- 

vl 

volutional, it can be conveniently viewed in this way. The above is one of 
two special cases of the fan-beam problem that had previously been solved [ 4 , 5 ]. 
An attempt was made here to use similar notation to simplify comparison. 

We ought to specify how g.(x') is found from g. ... As in the previous 

J ■ J 

paper [1], we approximate g.(x') by interpolation. If we sample N rays uni- 

J 

formly along a segment of length L of the line at right angles to the source-to- 
origin line, then 5x - L/(N - 1), and the v ray corresponds to 

X. = -1/2 + i 6X (57) 


Consequently, g.(x') is found by interpolation from g.,. and g,., , ... where 

J 1 j 0 + 1)J 


i' = L(x' + L/2)/6XJ 


(58) 
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In practice, of course, detectors would not be arranged on a line passing 
through the scanned space. The geometric transformations from a more distant, 
linear detector array to positions on the line passing through the origin are 
fortunately trivial (see figure 3). Such an array of equally wide detectors 
positioned behind the object being scanned would have to move in synchrony 
with the source, so as to always remain perpendicular to the source-to-origin 
line. There is great interest in scanning schemes which can instead use a 
fixed array of detectors. One such arrangement will be discussed in the 
next section. 






EVEN SAMPLING OF RAYS IN FAN ANGLE. 

Even sampling of rays in fan angle a can be achieved easily using a set 
of equally wide detectors arrayed on a sector of a circle with center at the 
source position. Curiously equal spacing of samples in fan angle can also be 
achieved when these detectors are instead placed on a circle passing through 
the source, with center at the origin (see figure 7). This follows from the 
fact that the angle at the center is just twice the angle at the source, and 
so equal angular spacing of detectors when viewed from the origin corresponds 
to equal angular spacing of detectors when viewed from the source. Such an 
arrangement of detectors has an advantage in that the detectors could remain 
stationary during scanning if the potential mechanical conflict between source 
and detectors could be solved. In any case, it is natural here to let 


£ = ot 

(59) 

So, 


D COS ajl = D COS a 

(60) 

Proceeding as in the previous section, we obtain (from equations 35 

and 36 ), 

2 it 

f(r, 4 >) = ^ f g"(a\e)0/K 2 ) de 

(61) 

, . +ir/2 

g" (a 1 ,B) = 1m J F [sin(a - a')] D COS a p(a,g) da 

e - 0 -w/2 

(62) 

where (by equation 25), 
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K 2 = r 2 + D 2 + 2rD sin (3 


(63) 


Next, we obtain the discrete approximation. 


= 4^r 2 9j(“') r 2 + D ' 2~ + ' 


2rD sin (g. - $") 66 j 

J 


(64) 


g., • 

a i J 


2 F., . D COS a- P-j 6a 

i l-l l 


(65) 


where 


a' - tan 


■[- 


r cos ( 3 . - <{)) 

--.J..- 

r sin ( 3 . - (|> 
j 


( 66 ) 


Once again p... is simply the i^* 1 ray sum in the j fan, while 6a is the (fixed) 

angular interval between rays in the fan. The angular interval associated 

with a particular fan, 6B-, is as defined before (equation 40). The filter 

3 

factors are: 


F k ” si 


k f 0 


(67) 


F o = 


-2 F k 

k 4 0 


( 68 ) 


Finally, one needs to detail the interpolation procedure for finding gj(a') from 

the discrete set of values g.,.. If N rays are sampled uniformly along an arc 

* 0 

XL 

of angle A, then 6a = A/(N - 1). The i ray then corresponds to 




(69) 


a. - -A/2 + i 6a 

l 

Consequently, 9j( a ') 1S found by interpolation from g^ tj . and g^., + jjj where 

i' = L(a 1 + A/2)/SuJ (70) 

This reconstruction method may be viewed as a pre-multiplication of the ray 
sums by D cos ct., followed by convolution, with a final post-multiplication by 
l/[r 2 + D 2 + 2rD sin(e. - <}>)]. This is the second special case of the fan beam 
reconstruction problem which had been solved previously [5, 63* 

Other fan-beam scanning geometries do not lead to such special case solu¬ 
tions however. Usually, a general linear operation is required. Fortunately, 
the method presented earlier allows one to treat arbitrary fan beam scanning 
geometries. Vie will study one in detail as an illustration. 
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A METHOD WITH UNIFORM SAMPLING DENSITY. 

Both of the scanning schemes discussed so far sample areas near the origin 
less densely than they do areas near the edge of the region of reconstruction. 
This can be seen when it is remembered that the ray sampling density is the in¬ 
verse of the Jacobian J [1] and that for fan beam scanning (equation 24), 

(71) 


(72) 


(73) 

The result of this variation in sampling density is that reconstructions will 
have somewhat better resolution (particularly in the radial, as opposed to 
tangential, direction) in outlying regions. While this effect is not very pro¬ 
nounced for fans that are fairly narrow, it is still of interest to investigate 
schemes providing uniform sampling density. That is, 

D cos a(|f-) = 1 (74) 

If this equation is integrated one finds. 


1 — 93 n 9ct 

0 - a* D cos 0 H 

Now for the first method (equation 38) 

i 

D cos a (||-) = COS 3 a 
While for the second method (equation 66) 

D COS a (|f-) = D COS a 
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* 



C 


z = D sin a 



This means that rays are spaced evenly in their perpendicular distance z from 
the origin (see figure 8). No convenient arrangenent of equally wide detectors 
will provide for sampling of the fan in this fashion, but clearly detectors of 
varying width arranged on either a straight or curved line can be used. Their 
width will increase with distance from the central detector. 

Now note that 


sin (a - a') = sin a cos a' - sin a* cos 


a 


(76) 


sin (a - a') = cos a‘(tan a - tan a') COS a 


(77) 


• / i \ / Q 2 - £ 

sin (a - a ) = -g—* 


Ty 


1 




2 - /qt-tj nr 


D 


(78) 


where $' = D sin a'. Proceeding as before (using equations 35 and 36), we get: 


2tt 


f ( r »*) = 4pr/ 9 " ( *’ 3) [D + r sTnnj^TTP 


dB 


(71) 


lim + D 

g-U'.B) ■ . J F ( 


e - 0' 11 /D 2 .-' ~l z W - (D 2 - A2) 


4 ‘ ) —J- rP (t,B)dA (80) 


The following identity v/as used for the outer integral, 

1 D 4 _ D 4 D 2 

* dT— .T = K^D^cos 2 a 1 = t'D +"r"sTn"TB - V)T 


( 81 ) 





Note that here V is the value of % for the ray inclined a' to the source- 
to-origin line (see figure 8), 


V = D sin a' 


(82) 


(This differs from the parallel ray scanning schemes presented in the previous 
paper [1]). 

Once again, a discrete approximation is required. 


f ( r » + ) : 4 ^ 2 [D + r s'in 


63 d 


(83) 


F £ I 1 


i J i vi 


p . . 6£ 


(84) 


where 


£* = D sin a' = 


Dr cos(g. - <|>) 

__—J- — «-«,■-— 

/r 2 + D 2 + 2rD sinfg. -J) 

J 


(85) 


The filter factors are 


F.,. = 


r-s- 

L^t: 


vzrr 

^rrj 


( 86 ) 


F . , = - 2 F.,. 

1 1 i i i' 11 


( 87 ) 


for i i i 1 




In this case, then, as for most scanning schemes, a general linear operator 

rather than a modified convolution must be used. The methods presented here 

permit the derivation of algorithms to deal with these problems. Note by the 

way that here the factors of sin (a - a') were split up in a similar fashion 

to how this had been done for the first two examples. This is not strictly 

necessary, since all three components can be accomodated as part of the filter 

function F (£,£') or F.,. if so desired. 

£ 11 





ANOTHER METHOD. 

The existence of an elegant method for reconstruction from parallel-beam 
data (equations 7 & 8 or 9 & 10) which uses derivatives of projection data and 
does not depend on arbitrary filter coefficients leads one to search for a 
similar expression for fan-beam reconstruction. Starting from the general 
form (equation 4) does not seem to lead to such a result. Instead one may apply 
partial integration to the form of the inner integral shown in equation 36, 

ir/2 

g"(a' ,n) = ]i n f F [sin(<x - a')] p(a,e) D cos a da (88) 

e+0 ' -ir/2 £ 


If one lets sin(s) = e, then 


a 1 -6 


g"(a',n) = [ / -—- p(a,e) D cos a da + 

0 ^ >2 sin 2 (a - a') 


a'+6 

~Y I p(a,3) D COS a da + 

a 1 -6 

i + " /2 i 

I-L- p(a,3) D COSa da] 

o , + slnMa - «■) 


( 89 ) 
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Or, 


a ‘ -6 


g"(a',n) = III [ ( ' - ■ 

-tt/2 sin2 (a - a') 


p(a,8) D - q ° S - --dal + 

COS(a - a 1 ) J 


(90) 


Y p(a 1 ,g)DC0Sa' + 


limf / +ir/2 

a +6 


a - a 


Sin 2 (a - a') 


P(a,8) 


Dcos a 

COS(a - a') 


That is, 


k 


_ n /„ p \ cos a 

^ COS(a - a' 


m(a - a 


a 1 -6 


•7r/2 


a '-6 


-tt/2 


_D_ r 3_ 

a - a' ) 8a 


[P(a»8) 


COS a 
COS(a - a 


da + 


p(a' ,8) COS a' + 


V-VT p( “' !l tszu 


COS a 

a - a 1 ) J 3 1 


+tt/2 


+tt/2 


a 1 +6 


a - a 


p(a,e) 


COS a 


(91) 


pi: 
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As 6 > 0 this becomes simply. 


+tt/2 


f D r / „\ COS a i , 

I sin(a - a') 3a a ’ COS (a - a') 

-tt/2 


(92) 


To summarize. 


f(r,<j>) = 


1 

4u 2 


I 


2tt 


g(a 1 ,e) 


D 


K 2 (r,<)>,&) 


de 


(93) 


g(a' ,b) 


- / 


7t/2 


-tt/2 


.- r -- 1 -[p(a,p) V u:> u tv ] da 

Sin(a - a j 3a LHV 'COS(a - a ) J 


COS a 


(94) 


Where, as before, 


K 2 = r 2 + D 2 + 2r D sin(6 - <J>) 


(95) 


, _ r cos(b - <t>) 

tana ~ [D + rsin(g - OT 


(96) 


The discrete approximation is 


f(r,<>) z — T g,-(a') 


D 


— 6B . 

4tt 2 j J ' K 2 (r,(j> ,3^) J 


j 


^ i cos (at) 

^ Sin (a. - a'., E p(d i’ B) COS(at -at,) 

j 1 1 1 1 


C°s(at ,) 

P( “i-1’ 6) a’,) 1 
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X L. 

Here corresponds to the left edge of the i tn detector, while marks its 

x U 

right edge (see figure 9). The ray-sum seen by the i Ln detector is p.. and its 

* J 

center is at a!. 
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CONCLUSION AND SUMMARY. 

The formulas for reconstruction from ray sums obtained by arbitrary 
sampling schemes were specialized to a system utilizing fan beams originating 
from sources on the circumference of a circle. It was found that one need 
not calculate the contribution of each ray sum to each point explicitly, but 
that the calculation does involve the application of a general linear operator. 

In special cases, this linear operator becomes space invariant by a manipulation 
of the integrals, and the superposition integral simplifies into a convolution. 

Two examples of this were shown, both corresponding to previously known solu¬ 
tions to the fan beam reconstruction problem for particular ray collection geometries 

To illustrate the utility of the new method, however, a third case was 
considered where the simplification does not occur. Previous techniques for 
finding reconstruction methods based on Fourier transforms cannot deal with 
it. Details of an algorithm were developed. The utility of the new methods 
for finding algorithms for arbitrary fan beam scanning schemes is therefore 
apparent. The introduction of uniform scanning coordinates in particular is of 
great importance in finding reconstruction methods for the general case. 
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